Phase transitions and crossovers in reaction-diffusion models with catalyst 

deactivation 
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The activity of catalytic materials is reduced during operation by several mechanisms, one of them 
being poisoning of catalytic sites by chemisorbed impurities or products. Here we study the effects 
of poisoning in two reaction-diffusion models in one-dimensional lattices with randomly distributed 
catalytic sites. Unimolecular and bimolecular single-species reactions are considered, without re- 
actant input during the operation. The models show transitions between a phase with continuous 
decay of reactant concentration and a phase with asymptotic non-zero reactant concentration and 
complete poisoning of the catalyst. The transition boundary depends on the initial reactant and 
catalyst concentrations and on the poisoning probability. The critical system behaves as in the 
two-species annihilation reaction, with reactant concentration decaying as t"^''* and the catalytic 
sites playing the role of the second species. In the unimolecular reaction, a significant crossover to 
the asymptotic scaling is observed even when one of those parameters is 10% far from criticality. 
Consequently, an effective power-law decay of concentration may persist up to long times and lead 
to an apparent change in the reaction kinetics. In the bimolecular single-species reaction, the critical 
scaling is followed by a two-dimensional rapid decay, thus two crossovers are found. 



I. INTRODUCTION 

Simple models for reactions of diffusing species have 
been intensively studied in the last decades [J 0, Q j such 
as trapping reactions, annihilation or coagulation of a sin- 
gle species and two-species annihilation. In low dimen- 
sions, these apparently simple models yield a wide range 
of nontrivial kinetic behaviors because mean-field the- 
ories (laws of mass action) fail. One-dimensional media 
may be a realistic description of the structure of catalysts 
with long and narrow pores, such as zeolites and porous 
oxides, or may represent step edges of two-dimensional 
surfaces where reactants preferrably adsorb. The former 
materials have inhomogeneous distributions of catalytic 
centers where the reactions take place, which motivated 
recent theoretical models. If the reaction is unimolecular 
with volatile products, then that inhomogeneity can be 
represented by trapping models 0, d, 0] , with the traps 
playing the role of catalytic sites. For other reaction 
mechanisms, the disorder in the spatial distribution of 
the catalytic centers was also considered in recent works 
0, H, [9I llOl [Tl| . In some cases, it leads to nontrivial 
dependence of the reaction rates on diffusion coefficients 
and catalyst density. Effects of catalyst geometry and 
conditions to improve the efficiency of a reaction process 
were also discussed in recent works [l3, [ij, [l3| . 

A phenomenon of great economic impact on industrial 
processes is catalyst deactivation, which is the reduction 
of catalytic activity due to blocking of active sites, dis- 
tortion or blockage of the porous structure, sintering of 
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metal particles etc [T^ [H, [T3] ■ One of the possible de- 
activation mechanisms is poisoning. It occurs when a 
reactant, product or impurity is strongly chemisorbed 
in the active sites, preventing further reactions but not 
affecting the diffusion of reactants and products along 
the pores. For instance, poisoning of metals by sulphur 
species {H2S, SO2, SOi'~ etc) is a problem in many 
catalytic processes, such as hydrogenation, methanation, 
Fischer Tropsch synthesis, steam reforming and fuel cell 
power production. In other cases, large aggregates are 
formed inside the catalyst pores; this is called coking or 
fouling and eventually leads to pore blockage. Due to its 
technological relevance, deactivation mechanisms are fre- 
quently included in models of catalytic processes [l^, [l^ ■ 
The main interest of those models is to predict the time 
evolution of the catalytic activity and the turnover fre- 
quency with a continuous flux of reactants, thus they are 
usually designed for a particular application. 

However, there is no systematic study of the effects of 
deactivation on the kinetics of simple reaction-diffusion 
models (e. g. trapping or annihilation). The aim of this 
work is to fill this gap by introducing one-dimensional 
models for catalyzed reactions of diffusing reactants with 
poisoning of catalytic sites. We will consider initial con- 
ditions with uniform distributions of reactants and no 
external flux. These conditions are far from those en- 
countered in real catalytic processes in porous media, but 
they are interesting as a first step to understand possible 
changes in those systems kinetics. 

We will study these models with scaling concepts sup- 
ported by numerical simulations. Unimolecular and bi- 
molecular (same species) reactions occuring upon contact 
of the reactants with the catalytic sites are considered. 
Poisoning is represented by the permanent blockage of 
those sites, with a given probability, immediately after 
the reaction takes place. Two limiting cases of this type 
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of model anticipate the presence of kinetic transitions: 
for very low deactivation rates and low initial reactant 
concentration, poisoning is negligible, thus that concen- 
tration decays similarly to the case without deactivation; 
however, for high deactivation rates and high initial reac- 
tant concentration, the catalyst will be rapidly poisoned 
and part of the reactants will not be consumed. This 
analysis raises the question of how these systems behave 
when deactivation rates, initial reactant concentrations 
and catalyst loadings change, and which types of transi- 
tion occur. 

We will show that two phases exist in those systems. 
In one phase, the reactant concentration continuously de- 
cay to zero, with the same asymptotic scaling of the main 
reaction. In the other phase, there is a non-zero reactant 
concentration at long times. In the boundary between 
those phases, the system behaves as in the two-species 
annihilation reaction, with reactant concentration decay- 
ing as t~^^^ in one dimension - the catalytic sites play 
the role of the second species in this case. Near a tran- 
sition point, we will show a long crossover to the asymp- 
totic scaling of the phase of decaying reactant concentra- 
tion. Thus, a correct interpretation of the main reaction 
kinetics becomes difficult, particularly with short time 
simulation data. For bimolecular reactions, the deacti- 
vation also enlarges the time interval of an intermediate 
two-dimensional scaling, where the concentration decay 
is faster than both the short-time and the long-time ones. 
These crossovers may also have important consequences 
in the interpretation of experimental results. 

The rest of this paper is organized as follows. In Sec. 
II we consider unimolecular reactions with inert products 
taking place at random catalytic sites (traps) subject to 
deactivation. In Sec. Ill, bimolecular single-species re- 
actions are considered and in Sec. IV we summarize our 
results and present our conclusions. 



which immediately occurs when A occupies the same site 
of C. With probability p, the catalytic site is deactivated, 
i. e. C is converted into an empty site (consequently, it 
does not affect diffusion of the other reactants). The 
reaction with subsequent poisoning corresponds to the 
reaction scheme 



A + C ^0, 



(2) 



where represents an empty site. Thus, the model may 
be viewed as a competition between the reaction ([T]), with 
probability 1 — p, and the reaction with probability 
p. This is illustrated in Fig. la. 
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FIG. 1: Schemes of (a) unimolecular and (b) bimolecular 
single-species reactions in a line with catalytic (filled squares) 
and non-catalytic (open squares) sites. In both cases, the 
step of a reactant towards a catalytic site is shown in the left 
and the possible outcomes are shown in the right, with the 
associated probabilities. 



II. TRAPPING REACTIONS WITH TRAP 
DEACTIVATION 

In this model, defined in a one-dimensional lattice, we 
consider unimolecular reactions with volatile products, 
which occur when the reactant is in contact with the 
catalyst. Initially, catalytic sites (C) are randomly dis- 
tributed with coverage ctq and are immobile. The other 
sites are labeled as non-catalytic or empty sites. At 
t = 0, a concentration of reactants (species A) is ran- 
domly distributed through the non-catalytic sites - con- 
sequently, the condition (Tq-I-^o < 1 is obeyed. For t > 0, 
reactants A diffuse with coefficient D, i. e. each reactant 
executes an average of D random steps to neighboring 
sites per unit time (the lattice parameter is the length 
unit), and they obey the excluded volume condition. 

The unimolecular reaction corresponds to a trapping 
reaction [l[ in the form 

A + C^C, (1) 



As time evolves, the catalyst coverage and the reactant 
concentration are respectively denoted by a (t) and 9 (t) . 

We performed simulations of this model in lattices of 
length N = 2^° (more than one million sites), with sev- 
eral different values of the parameters ao, Oq and p. For 
each set of parameters, typically 10'^ different realizations 
were averaged. Diffusion coefficients of species A were 
_D = 1 in most cases. Maximum simulation times were 
up to 10^ units under these conditions. 

A. Long time scaling 

For p — 0, we recover the well-known trapping problem 
with static traps. Since there is no deactivation, all the 
reactants initially placed in a segment between two con- 
secutive traps are annihilated in one of the edge traps. 
For a fixed segment size L, the decay is a simple ex- 
ponential in t^/^/L. However, the average over (Poisson 
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distributed) segment lengths is dominated by trapping in 
large segments which leads to the Donsker-Varadan 
11 result 
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(3) 



In the 



where p is the trap density and a 
case p = 0, we have p — cto- 

For p small but non-zero, the concentration of reac- 
tants rapidly decays to very low values, while only a 
fraction of the initial catalytic sites is deactivated. The 
annihilation of a concentration of species A leads to 
poisoning of a density p^o of catalytic sites. Thus, at long 
times, the trap density is cr = do ^ P^o and the reactant 
concentration 9 is much smaller. In this case, reaction 
P]) is asymptotically dominant, thus a decay as in Eq. 
([3]) is expected. 

Simulations with small values oi p confirm this result, 
as shown in Fig. 2a. They also show another interesting 
and somehow surprising feature: while results for p = 
show significant corrections to the dominant scaling 
(Eq. [3|), the results for non-zero p show that stretched 
exponential decay with negligible corrections (the data 
for p = 0.1 in Fig. 2a fits a straight line in three time 
decades). The corrections for p = were observed in 
previous simulation work 'F], and for this reason it took 
many years for the original analytical prediction in two 
and three dimensions Q to be confirmed numerically Q . 

From the above discussion, we expect the same asymp- 
totic behavior for any initial condition where ctq > P^o- 
On the other hand, if uo < P^Oi a- finite reactant concen- 
tration 9 = Oq — gq/p is found at long times, while the 
trap density continuously decays. This regime is dom- 
inated by reaction with excess concentration of A 
relatively to C. Fig. 2b confirms this result for = 0.4, 
ctq — 0.15 and p — 0.7, which gives a remaining density 
of reactants 9 « 0.186. 

A transition is found when 
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which corresponds to a surface in the (ctqi^OiP) space. 
In this case, the reactant concentration and the catalyst 
coverage decrease by the same factors at all times. 

The continuous increase of the distance between con- 
secutive catalytic sites certainly rules out the arguments 
that lead to the result in Eq. Instead, in a transi- 
tion point this system resembles the two species annihi- 
lation reaction of Eq. ([2]) with balanced concentrations 
of species A and C. According to previous work P, [3, 
in the reaction ^ with equal initial concentrations po of 
A and C in one dimension, the reactant concentration p 
decays as [l|, [loj 



1/2 



[iD,-,Dc)tt'^' 

where Z)^ and Dc are diffusion coefhcients. Note that 
this scaling law does not change if one species is static, 
which is the case of the catalyst C in our model. 
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FIG. 2: Reactant concentration as a function of time for the 
unimolecular reaction with ctq — 0.15, Oq = 0.4: (a) Exponen- 
tial decay for p = (solid line) and p = 0.1 (dashed line), (b) 
Decay towards a non-zero value for p — 0.7. 



The decay as is confirmed for several transition 

points, as shown in Fig. 3. However, the dependence 
with the initial reactant concentration is different from 
the square-root law of Eq. ([5]) : for fixed a^, 6 is propor- 
tional to 6*0. This occurs because the kinetics leading to 
Eq. (I5|) is a balanced annihilation of AC pairs, but in our 
model the annihilation of an AC pair requires trapping 
of 1/p reactants (A) in the average. Thus, po must be 
replaced by 9o/p = ^o^/fo in Eq. ([5]). 



B. Crossover scaling 

Here we focus on the case where the system is near the 
critical surface (Eq. [4]) but p < pc = (Jo/9o, i. c. there 
is a continuous decay of the reactant concentration. At 
short times, this decay follows the critical power law of 
Eq. (O, but it crosses over to the stretched exponential 
behavior (Eq. [3]) at long times. This is illustrated in 
Figs. 4a and 4b for fixed values of ctq and 6*0 and various 
P- 

One interesting point is that, even with one parameter 
being 10% distant from the critical value (e. g- p = 0.34 
in Fig. 4a), the concentration decay for 9 > 10""* 
still resembles a power law. The exponent of this ap- 
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FIG. 3: Critical decay of reactant concentration in the uni- 
molecular reaction model with ao = 0.15: 6o — 0.75, p — 0.2 
(circles); do = 0.5, p — 0.3 (squares); = 0.3, p — 0.5 (dia- 
monds); 6o — 0.1875, p = 0.8 (triangles). The solid line shows 
a decay as t~^^^. 
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FIG. 4: (a) Short time and (b) long time decays of reactant 
concentration in the unimolecular reaction model near the 
critical boundary with ao = 0.15, ~ 0.4: p — 0.34 (circles), 
p — 0.35 (squares), p = 0.36 (diamonds), p — 0.37 (triangles); 
critical point is at p = 0.375. In (a), the dashed lines have 
slopes —1/4 and —1. 



10^ <t< 10^, as illustrated in Fig. 4a. Closer to tc, ap- 
parent power laws are found for longer times. This is a 
surprisingly long crossover due to the remarkable differ- 
ence between the scaling of Eq. ([5]) and the asymptotic 
one of Eq. 

At long times, the remaining density of catalytic sites 
is proportional to pc — p. Eq. ([3]) gives a reactant con- 
centration 



^ exp \-biDt)'/''] ,b^{pc- pf 



(6) 



However, the scaling oi 9 on pc ~ p is difficult to be con- 
firmed numerically (e. g. in Fig. 4b) because the data 
close to Pc shows huge corrections to the (noncritical) 
asymptotic decay of Eq. [31 

From the above result, we can estimate the time in 
which the apparent power law decay crosses over to the 
exponential one. Matching the time scaling in Eqs. ([6]) 
and ([5]) (the latter with no amplitude depending on pc — 
p), we obtain a crossover time of order 



(7) 



This helps to understand the long crossover in this sys- 
tem: for the values of p shown in Figs. 4a and 4b, we 
expect tc to range approximately from 10^ to 10^. Ex- 
perimentally, it has an important consequence: the re- 
action kinetics may be incorrectly identified if one ob- 
serves the time evolution of the concentration decay (or 
the turnover rate) in a restricted time range and disre- 
gards the effects of deactivation. 



III. ANNIHILATION REACTIONS WITH TRAP 
DEACTIVATION 

The second reaction mechanism studied in this work 
is bimolecular with a single species. The model corre- 
sponds to annihilation reactions of this species limited to 
catalytic sites of a lattice. 

Again we consider initial random distributions of static 
catalytic sites C, with coverage (Tq, and diffusing reac- 
tants A with initial concentration Oq. The bimolecular 
reactions between two A which meet at a catalytic site is 
represented by 



A + A + C^C, 



(8) 



It corresponds to the well known one-species annihila- 
tion model {A + A ^ 0) restricted to the set of catalytic 
sites, as in Ref. 0. The poisoning is also represented 
by the permanent blockage of a catalytic site after a re- 
action, and occurs with probability q. It corresponds to 
the reaction 



parent power law depends on the chosen time range, 
varying from —1/4 for 1 < t < 10^ to almost —1 for 



A + A + C-^O. (9) 
The competitive model is illustrated in Fig. lb. 
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A. Long time scaling 



When 9 = 0, we have the model proposed by Oshanin 
and Blumen [7| , in which the density of reactants decays 



l/iDt) 



-1/2 



(10) 



in one dimension. This dominant contribution to the de- 
cay does not involve the catalyst coverage, which can 
be explained as follows. The characteristic time for 
two consecutive reactants to meet at the same site is 
taa ~ 1/ {DO"^). After they meet, the closest trap 
is inside a neighborhood of size I/cto. Thus, the typ- 
ical time for those reactants to meet at that trap is 
taac ~ 1/ (-Dctq^)- At long times, 9 <C ao, thus we have 
TAAC ^ TAA and the long-time scaling is dominated by 
TAA, which leads to Eq. (fTU]) . However, the catalyst cov- 
erage plays an important role at intermediate times, as 
will be shown in Sec. IIII Bl 

For small g, we expect the same scaling as in Eq. ([T0)|. 
At long times, 9 ^ 9o, thus the density of annihilated A 
particles tends to and a density q9o/2 of catalytic sites 
will be deactivated. Thus, the final catalytic coverage is 
(To — q9o/2. Simulation results shown in Fig. 5 confirm 
these predictions. 
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FIG. 5: Decay of reactant concentration far from the critical 
boundary for the bimolecular single-species model with ctq = 
0.1, do = 0.3: q = (circles) and q — 0.1 (squares). The 
dashed line shows a decay as t~^^'^. 



We expect a transition between a phase of decaying 
reactant concentration and a phase of non-zero long time 
concentration when 



concentration at the transition point is that of a three- 
particle annihilation model. Since the catalyst coverage 
does not play a role in the decay for q = (Eq. [TOl), we 
might think that the time taa governs the scaling at the 
transition. However, this is not the truth. Similarly to 
the unimolecular case of Sec. II, there is accummulation 
of C at some regions where the initial concentration of A 
was locally small. Also, there is accummulation of A at 
neighboring regions, where its initial concentration was 
locally high. Thus, the coarsening is again dominated 
by the annihilation time of neighboring domains, which 
scale as 1/9'^ [l|, [l^, and not by the time taa- Thus, 
this transition also scales as in two-species annihilation 
model, although three particles are involved in the reac- 
tions IHlandin]). 

These conclusions are supported by simulation results. 
In Fig. 6a, we show the average number of consecutive 
A, (Sa), and of consecutive C, (Sc), as function of time, 
in one point of the critical surface. They scale approxi- 
mately as 



{Sa) ^ {Sc) ^ t 



1/4 



(12) 



Fig. 6b shows the decay of the reactant concentration 
approximately as (Eq. ^ in three points of the 

critical surface. All these results are typical of the two- 
species annihilation problem [i| • Note that the total size 
of domain A (C) includes the reactants (catalysts) and 
the empty sites between them, thus this size increases as 
t^/^ (in the literature, this total size is the one usually 
referred to as the domain size). 

In Fig. 6b, we observe that the critical 9 depends 
on the initial concentration in the same way as the uni- 
molecular reaction, i. e., with p in Eq. ([5]) replaced by 
9olq - Oo^/ao. 

For CTo < the catalyst coverage vanishes asymp- 

totically, while a finite number of reactants is present. 
The final (non-consumed) reactant concentration is — 
2(Ja/q. These results were also confirmed by simulation. 



B. Crossover scaling 

The above analysis suggests that a crossover between 
two power law decays (i~^/'* and f~^/^) would be ob- 
served near a transition point. However, again the situ- 
ation is much more complex. Ref. Q analyzed the prob- 
lem with (7 = (no deactivation) and showed that, for 
small catalyst concentration, there is a long time interval 
with effective two-dimensional decay of the concentra- 
tion: 



ao = q9a/2. (11) 

In this critical surface, reactant concentration and cat- 
alyst coverage decrease by the same factors as time in- 
creases. 

The three-particle reaction ^ and an analogy with 
the case of Sec. II could suggest that the scaling of the 



e-lnt/t. (13) 

Small catalyst concentration is actually the case as we 
approach the critical point. Consequently, not only two 
but three scaling regions are expected: for short times, 
the critical decay as while the orders of magnitude 

of (7 and q9/2 are not very different (see Eq. [TT|) : for 
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FIG. 6: (a) Average number of consecutive reactants (circles) 
and consecutive catalytic sites (squares) in a critical point 
of the bimolecular single-species reaction model; cto = 0.4, 
9o = 10/11, q = 0.88. The dashed line shows a scaling as t^^^. 
(b) Critical decay of the reactant concentration in the same 
model with ao = 0.1: q = 0.4, 9o — 0.5 (circles); q = 0.5, 
00 = 0.4 (squares); q = 0.8, 9o = 0.25 (diamonds). The 
dashed line shows a decay as t'^^"^. 




FIG. 7: Reactant concentration near the critical point of 
the bimolecular single-species reaction model with ctq = 0.1, 
9o = 0.8: q = 0.21 (circles), q = 0.22 (squares), q = 0.23 
(diamonds), q — 0.24 (triangles); critical point is at g = 0.25. 
Dashed lines indicate three types of power-law decay present 
in this system. 



intermediate times, the much faster power law decay of 
Eq. [T31 for t ^ l/tToo^ [3], where cToo is the asymptotic 
concentration of the catalyst, the asymptotic scaUng (fTO|) 
is found. Since l/coo^ ^ {Qc ^ ^ near criticality, the 
deactivation significantly enlarges the time window of the 
intermediate scaling. 

These results are illustrated in Fig. 7 for several val- 
ues of q. For q — 0.21, (16% far from criticality), a very 
long intermediate regime is observed, and the asymptotic 
decay of Eq. PH]) is found only for t ^ 10^ and concen- 
trations below 10"*. Note that this decay is not depen- 
dent on initial concentrations, as shown by the data for 
q = 0.21 and q — 0.22. Closer to the critical point, the 
asymptotic decay is not observed even with concentra- 
tions below 10~^. Compared to the unimolecular reac- 
tions (Sec. II), finding the true reaction kinetics from 
reactant or product concentration data is much trickier 
in this case due to the presence of two crossovers. 



IV. DISCUSSION AND CONCLUSION 

We studied reaction-diffusion models in heterogeneous 
one-dimensional lattices with catalytic and non-catalytic 
sites and possible deactivation of catalytic sites by poi- 
soning. These models show a transition between a phase 
with continuous decay of reactant concentration and sur- 
vival of catalytic sites, and a phase with asymptotic non- 
zero reactant concentration and complete poisoning of 
the catalyst. The critical systems behave as the two- 
species annihilation model with stoichiometric concen- 
trations of reactants. For each reaction mechanism, the 
transition boundary depends on the initial concentrations 
of reactant and catalyst and the probability of deactiva- 
tion. We analyzed models of unimolecular and bimolec- 
ular single-species reactions, and showed that, in a tran- 
sition point, the reactant concentration decays as in the 
two-species annihilation reaction. In the decaying phase 
but near a critical point of the unimolecular reaction, a 
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long crossover to the asymptotic scaling is found. Thus, 
short time simulations show effective power-law decays 
instead of the stretched exponential, which suggests a 
change in the reaction kinetics. The situation is much 
more complex with the bimolecular single-species reac- 
tion because the small catalyst density near a critical 
point leads to an apparent two-dimensional decay of re- 
actant concentration. Thus, two crossovers are observed 
and the asymptotic decay is found only at extremely long 
times, even with one parameter being more than 10% far 
from criticality. 

Our results resemble transitions in other competitive 
reaction-diffusion models which were not explicitly inter- 
preted as effects of catalyst deactivation. Sanchez and co- 
workers studied the case where reaction A + B ^ B com- 
petes with B + C ^ C (double trapping) oi B + C ^ 
(trapping plus annihilation) [28l.[29j. showing that the ini- 
tial concentrations and the reaction rates also determine 
the asymptotic behavior. Further work on these and 
related reactions also considered anisotropic or anoma- 
lous diffusion of reactants, different forms of reactant 
input, and power-law decays of reaction rates, which 
lead to a wide range of interesting kinetic behaviors 



[13, HTUp, m, m, m . For instance, the model of Yuste 
et al [35j| for trapping by subdiffusive traps with vanish- 
ing density shows a transition with scaling similar to our 
unimolecular model. 

Despite the simplicity of the reaction-diffusion models 
presented here, experimental applications are possible. 
For instance, the depletion zones of a trapping model 
with a single trap were already observed in photobleach- 
ing of fluorescein dye by a focused laser beam [H, [sj , 
and models with competitive reactions were already used 
to represent patterns of a reaction between Cr^~^ and 
Xylenol Orange [3^]. Thus, we believe that our work may 
also have applications to systems with catalyzed reac- 
tions and deactivation by poisoning, particularly to help 
the interpretation of kinetic data. 
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